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Abstract 


An improved understanding of the large-length-scale capillary flows arising in a low-gravity 
environment is critical to that engineering community concerned with the design and analysis of 
spacecraft fluids management systems. Because a significant portion of liquid behavior in 
spacecraft is capillary dominated it is natural to consider designs that best exploit the 
spontaneous character of such flows. In the present work, a recently verified asymptotic analysis 
is extended to approximate spontaneous capillary flows in a large class of cylindrical containers 
of irregular polygonal section experiencing a step reduction in gravitational acceleration. Drop 
tower tests are conducted using partially-filled irregular triangular containers for comparison 
with the theoretical predictions. The degree to which the experimental data agree with the theory 
is a testament to the robustness of the basic analytical assumption of predominantly parallel flow. 
As a result, the closed form analytical expressions presented serve as simple, accurate tools for 
predicting bulk flow characteristics essential to practical low-g system design and analysis. 
Equations for predicting comer wetting rates, total container flow rates, and transient surfaces 
shapes are provided that are relevant also to terrestrial applications such as capillary flow in 
porous media. 


Introduction 

A significant portion of liquid behavior in spacecraft is capillary dominated: as for liquid 
propellants, cryogens, thermal fluids, and wastes. It is therefore natural to consider designs that 
best exploit the spontaneous character of such capillary surfaces and flows. Fundamental insights 
relevant to a variety of “characteristic” capillary flows continue to be sought that will further 
enable rapid and accurate predictions of important features of the fluid behavior. 

Recent investigations have successfully demonstrated asymptotic techniques for the solution of 
capillary flows in containers with interior comers (see Ref. 1 and references contained therein). 
Such flows are “characteristic” of low-g fluid behavior in propellant tanks employing vanes or 
baffles for passive fluid positioning. In several cases, key features of the flow such as volumetric 
flow rate and surface shape may be determined in closed form — valuable tools to the designer of 
space-based fluids management systems. 

The problem of sudden capillary rise (imbibition) in containers with interior comers is common 
to drop tower tests and serves as a model problem for liquid/tank filling, draining, response to 
thruster firing, docking, etc. In this paper an asymptotic analysis is briefly reviewed that has been 
successfully demonstrated to predict such flows in simple cylinders of rectangular and equilateral 
triangular section. The analysis is then generalized and extended to approximate flows in 
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cylindrical containers of irregular polygonal section with special attention paid to the 
development of a global similarity solution. Fully transient 3-D surfaces derived from the closed 
form expressions are presented for comparison with experiments. The results of simple drop 
tower tests employing irregular triangular containers are reported and the theoretical predictions 
are used to benchmark the salient features of the flow such as time dependent liquid column 
length and flow rate. Lastly, the value of the design relationships is briefly summarized. 

Review of Ref. 1 : Single Comer Solution 


Detailed comparisons between experiments and theory have demonstrated the validity of the 
assumption that the flow throughout the container is controlled by the local capillary flow in the 
comer 1 . Assuming a wetting fluid and locally parallel flow [(HAS) 2 « 1], the dimensionless 
leading order governing equations simplify to the nonlinear lubrication PDE 

h x - 2h: + hh„ ( 1 ) 

where //(-, r) is the dimensionless height of the meniscus measured along the bisector of the 
comer at location z and scaled time r(see Figure 1 for notation). This governing equation implies 
that the surface may be approximated as a construct of circular arcs in the cross-flow plane (x-y 
plane), and, once h(z, t) is determined, the entire 3-D transient surface may be determined from 

S(y,h)=h(l + f)-(f 2 h 2 - y 2 tan 2 a) !/2 (2) 

The parameter/ is the measure of interface curvature (driving force) and is given by 



where 0 and a are the contact angle and comer half angle, respectively. The static contact angle 
boundary condition is correct to leading order and is thus applied at the contact line. Dynamic 
contact angle effects can be relegated to higher order, as the predominant flow direction is 
parallel to the contact line. Eq. (1) is in a general class of PDEs that emerge in nonlinear heat 
conduction 2 and foam drainage 3 . The problem of sudden capillary rise (imbibition) applies 
constraints /i(0,t) = 1, /z(l,r) = 0, and conservation of mass to Eq. (1). Transforming Eq. (1) by 
h = F(i)) and 1 ] = z(2zy ]r2 yields the similarity equation 

FF rvi +2F 2 +l 1 F^0 (3) 

subject to F(r\, ip ) = 0 and F n (rj, ip ) = -1/2. Eq. (3) is invariant under the transformation F = A 2 F* 
and r)= X if , which may be exploited so that the meniscus tip may be conveniently located at if 
= 1, simplifying numerical calculation and producing the system 

F + F; v+ 2F^ + ^Ft=0 (4) 


subject to F f ( 1) = 0 and F f ^(l) = -1/2. The numerical solution of F*(rf ) is provided in Figure 2. 
Because h is scaled by H, which is known, A 2 = [F(0)]' I/2 , and with F*(0) = 0.345, A = 1.702. 
As mentioned above, the arbitrary and dimensional length L scales ~, and the dimensionless time 
r is given by 


T 


H <7 Esin 2 a 
— — i 1 

2L- n f 


(5) 


where t is dimensional time, a and p are the fluid surface tension and dynamic viscosity, 
respectively, and sin 2 a/ f is the geometric ratio of capillary driving force (1 If) to viscous 
resistance ( 1/ F t sin 2 Ot ). F, is a weak function of 0 and a (see Ref. 1 , Fig. 6 for exact value). 
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The solution of Eq. (4) for a single interior corner of infinite length provides important design 
quantities such as liquid column length L and flow rate Q as functions of time. These quantities 
are provided below in dimensional form: 

L = 1.702 G U2 H U2 t >/2 (6) 

Q = 0.349 / 2 F An G' r -H' n -r ul (7) 


with 


r cos# cos(a + 6) , 

F An = T— — + « + » 

sin a 2 

and where H is a known constant height condition at - = 0. G is given by 

oF sin 2 or 


G =- 


Hf 


( 8 ) 

(9) 


With eqs. (6) through (9) low gravity containers may be sized, optimized, fluids selected, or flow 
times predicted. Such quantities, which can be rapidly computed with a calculator, are accurate 
to ±6% for perfectly wetting fluids' and represent an improvement over previous design 
relationships that used comer friction factors and weighted capillary pressures 4 . Even the flow in 
the vicinity of the bulk receding meniscus may be determined. Despite the fact that in this region 
the fundamental assumptions of parallel flow and 1-D interface curvature no longer apply, 
knowledge of the fluid removal rate from the bulk, which is well-described by the individual 
comer flows, enables the prediction of the bulk meniscus recede rate. The general concept 
behind the “global flow” solution was introduced in Ref. 1 and reveals similarly that the bulk 
meniscus axial location (for cylinders of regular polygonal section) obeys the relation 

Z oc{ a tjn) 11 (10) 

the constant of proportionality being a function of specific container size and shape. 


Extension to n-Sided Polygonal Containers 


General and Global Flow Characteristics 


The accuracy of the closed form expressions of eqs. (6), (7), and (10) have been established 
experimentally in drop tower tests employing equilateral triangular cylinders and rectangular 
cylinders of various aspect ratio 1 with 0=0. With this success, the solution approach may be 
generalized and extended to a large class of irregular polygons without re-entrant comers 
provided the Concus-Finn condition 1 ' 5 , 9< n!2-a jy is satisfied. A generalized irregular polygonal 
cylindrical section is shown in Figure 3 (left) to provide the necessary notation for the results to 
follow. The solution is limited to containers where the liquid “rises” independently in each 
comer that satisfies the Concus-Finn condition. The solution does not apply to irregular polygons 
where a single interface may span two or more comers. An example of such a container is 
depicted in Figure 3 (right) which may be only approximately analyzed using the techniques 
presented herein. 

Employing the method of de Lazzer et al. 6 to evaluate the capillary pressure in the infinite 
container it is possible to show that 


Lj = 1.702 Gj ,2 Hj /z t U2 
Qj= 0.349// F M) G J U2 H/ ,2 r' n 


( 11 ) 

( 12 ) 
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where 


H 


Q,o ,=HQ, 

j = I 



aF tj sin 2 cr 


fj 


' cos 6 j 
since 

J / 


P n cos 9 
2 /; S 



44 X 1 

i/:~ 

1- 

1 — ^^7 



P~ cos“ 9 


_ 

\ n > 



( 13 ) 


A;? j 



2 =ix, 

y=> 

cos0cos(ce; +0) ;r 

^ +«, +0 

sin a ; 2 1 


where j denotes the j" 1 of n interior corners. P„ and A„ are the perimeter and cross-sectional area 
of the n-sided polygon, respectively. It is useful to note from eqs. (11) and (12), that for a given 
polygonal vessel and fluid, 

L ; °c ] - since, (14) 

Qj « (l-sinorpfcoto^. -n 1 2 + a y ) (15) 

To give an example of the utility of these expressions, for a 45-45-90 isosceles triangular 
section the ratio 


r 4 s r l -sin 22.5° 210 

Ig o 1- sin 45° 


(16) 


Therefore the liquid column in the 45° (a= 22.5°)interior comer is always about 2.1 times longer 
than the liquid column in the 90° comer. A similar comparison may be made for Q } . The purely 
geometric relationships in eqs. (14) and (15) may be used to simplify hand calculations; the 
calculation of L, and Q, is necessary in only one comer, the Eq. (14) and (15) relationships are 

then used to form ratios by which to quickly determine L, and Qj in the remaining comers. 


The results of eqs. (11) and (12) above are correct provided the coordinate origins for the 
individual comers are known. These are the locations of the constant pressure, constant height 
location (as determined by de Lazzer et al. 6 ) boundary condition for the individual comers. 
Assuming that the coordinate origins for each comer coincide, a global similarity solution for the 
flow throughout the container is possible that is also fortunately consistent with the initial 
condition of a typical drop tower test. The pertinent details of this analysis are provided below. 
Further discussion of the “common origin” assumption is provided in a subsequent section. 


In Figure 4 is sketched an irregular container section (left) with contrived interfaces shown in 
profile along the B-B plane of symmetry (right). The comer rise rates are different in comers of 
different Cf, and the assumed common origin is identified where the constant height condition 
hj(0, r) = 1 for each comer applies. Because the details of the comer flows are known precisely in 
the domain 0 <Zj< ( Zu p )j , the objective of this “global” analysis is to provide the interface shapes 
within Zb < Z< 0 and compute the bulk meniscus location Zb as well. The assumption of parallel 
flow in the individual comers is valid throughout the domain, z b < Zj < ( z,( P )/. note that {z,i P )k = L k . 
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The analysis begins with a dimensional mass balance about the origin, noting that the volume of 
fluid drawn up the comers (z > 0) is equal to the volume of fluid removed from the bulk (z < 0); 
the volumes of the cross-hatched regions above and below - = 0 in Figure 4 are equal. Therefore, 
noting Zb < 0, 


)1 '-IW 


'L}r v h:d; = -A.z b --ZjF tl h;d; 

; =! o M z b 

where F Aj = f/F Anj . Combining terms yields 

77 

4,*+X J^*;*=o 


(17) 


;=■ s. 


Nondimensionalizing h, and z by Hj and L, and introducing similarity functions hj = F/ rj/i, X] } = 
z(2xj) m , r),i P j- Ziipji 2tj)' ,c , and 77^(2 r,) 1/2 , as accomplished for the governing similarity ODE 
Eq. (3), Eq. (17) is transformed to 


2, F„Hl( T ' 


1/2 


AJ“J 


» 

7=1 f lbr 


J 

r 


Thipj 


i F i d i ,= o 

QbrlTr'Tj)' 1 * 


(18) 


The global similarity solution naturally assumes that the bulk meniscus location behaves 
according to the t m law and employs a reference time scale based on arbitrary reference angle a ,. , 
interface height H„ and interface curvatures /■ reminiscent of Eq. (5), namely 

H, a f),siirg, . 

r 2L 3 [I f. 

Thus, noting Hj = R/fj from Eq. (13), 


\ 

It 

1/2 

fj since. 

F T 

tr 

UJ 

f r sina ; 

KJ 


Note that the weak functional dependence of F,y is considered despite the fact that this parameter 
is bound 1/8 < < 1/6 for all possible values of 6 and a,. Substituting invariant transform 

functions r\ } = X, rtf and Fj = Xf F/, Eq. (18) transforms to 

A+tf-fi it 1f = <= 0 < l9 > 

j- { ” hr j ff+Bj 


where 


and 


w+ _ ntF'; 2 sina r 
" /, 


R = f> 

1 F, 1/2 sin« ; 


ij j 

At this point it is clear from Eq. (19) that the integrand is independent of comer j, thus Xj = X, = 
X and, after further rearrangement noting Hj = R/fj and F AJ =ff F AnJ , Eq. (19) simplifies to 


UbrA, ■ y 

R 2 


J (X 2 F + ) 2 drj + = 0 


( 20 ) 
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which must be solved implicitly for rj* r once f* is computed from the invariant similarity ODE, 
Eq. (4). As previously stated, the solution for F* shown in Figure 2 gives F*( 0) = 0.345... which 
yields A = 1.702... It may also be shown that as rf becomes increasingly negative, F*{ if ) 
asymptotically approaches -0.465... Thus, from Eq. (20) it may be shown that as cost? - sina ; 
approaches 0 for a given comer the volumetric flow along the comer contributes negligibly to the 
value of rf£. . 


Computed values for ff^ r from Eq. (20) are listed in Table 1 for a variety of familiar container 
cross sections. In the limit I Tj^ r B } I « 1 it may be shown that 



-0.4103^F Al!/ /v /2 (cos0 -sincr) 


7=1 



( 21 ) 


which is a leading order expression revealing most clearly the general dependence of on the 

geometry of the system. Values of Jfbr computed from Eq. (21) are also listed in Table 1 where it 
is seen that errors incurred are < 2% for most of the containers listed. Errors using Eq. (21) 
increase significantly however for small comers O' and slight curvature cos0 - sinc^, where 
« 1 is no longer true. 


With rj h + r known, the dimensional bulk meniscus location is found from 


^= 1 . 702 ^ 


/ \i/2 

f Ra ' 






( 22 ) 


and because t]£ r < 0, z h is of course increasingly negative in time. Cast in this form it is insightful 
to note that the liquid column length in each comer from Eq. (11) is 


Lj =1.702 F" 2 (cos 6- sin a y )| 


Ro 

i 


0/2 


) 


— K,. t 


(23) 


and the total active interface length for each comer may be determined by 


-Lj — z b =\ .702 


Ra 
1 

p 


i 


/j (Fj / \cosd-sma j )-Tf*)=K tolJ t ul (24) 


Thus, the global similarity solution is given by Eq. (23) for the comer flows, Eq. (22) for the 
receding bulk meniscus, and Eq. (24) for the total “hydrodynamically active length” along each 
comer. 


Discussion on Global Solutions 


Included in Table 1 is a term Q geom that represents the geometric dependence of the total flow rate 
along the interior comers of the container, 

Q„„ » 100 ^ ,! (cose-sina,) « Q,„ (25) 

7=1 


NAS A/CR— 2001-210900 


6 


Q geom is computed for each cross-section holding cross-sectional area fixed. It is important to note 
that Qgeom and thus Q tot tracks with Tf £. . It is also important to note that Q gcom changes 
dramatically with container shape for fixed cross-sectional area. For example, the equilateral 
rhombic sections computed in Table 1 reveal at least a 23-fold increase in volumetric flow rate as 
the acute angle is reduced. These trends illustrate that for fixed cross-sectional area, containers 
with a larger number of smaller comer angles will transport significantly larger amounts of liquid 
by capillary forces, i.e. a porous material with 30-150 ( = 30°- 150°) equilateral rhombic channels 
will transport 9-fold the liquid of a porous material with 90-90 equilateral rhombic channels with 

the same total cross sectional area! It can be shown that Q WJ ^ OC l/ 4 and thus approaches zero 

slowly for vanishingly small values of a as indicated by reduced value of Qgcom listed in Table 1 
for a 0.11-179.9 equilateral rhombus. Qgcom values for the regular n-gons approach zero with 
increasing n as expected. 

Detailed Interfacial Results 


Perhaps the most important design quantities summarizing such flows are provided above in eqs. 
(11), (12), and (22). However, the entire surface profile of the liquid throughout the container 
may be computed if desired. The solution follows from the generalized global similarity solution 
and is applicable at long times throughout the container, despite the fact that both the flow and 
interface shape are not known in the neighborhood of the bulk meniscus. By approximating the 
global similarity solution for the meniscus centerline height in each comer by the polynomial 

hj =// ; (1-0.57 \rf -0.429r7 +2 ) (26) 


where H ■ is given by Eq. (13) and 


7]* =0.587 


Pfj 


xl/2 


zt~ 


(27) 


H j(F-) j sin~ GCj J 

subject to the constraint ff br Bj <rf + <1, the 3-D transient interface in each comer may be 


computed via 

Sj=hj (\ + fj) + (h~ ff-y ]) 1 ' 2 


(28) 


where 

\yj\ ^ hj fjCO $( 0 Cj + 9) 

Samples of computed surfaces using this simple approach are presented in Figure 5 for 45-45-90 
and 30-60-90 triangular cylinders. Note that the column length for the 45° comer is about 
2.1 times that of the 90° comer as predicted by Eq. (16). These results extend previous solutions 
for equilateral polygons 7 . 


Impact of Approximations 

Essential to the analyses above is the assumption that the constant height condition for each 
comer occurs at a common c-coordinate origin. The validity of this assumption has not been fully 
established, but it will be demonstrated experimentally that the general solution approach is 
rather insensitive to the precise location of the coordinate origin(s). The common origin 
assumption permits a global similarity solution that perhaps best captures the initial condition of 
typical drop tower tests. However, there are certainly unphysical characteristics of the solution at 
small times. In Figure 6 is shown 3 frames of a drop tower test the details of which will be 
provided in the experiments section. For the similarity solution the constant height location Hj for 
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the subject comer flow occurs at what is assumed to be the "-coordinate origin (z = 0) for each 
comer flow. It is important to note that this location does not coincide with the initial interface 
origin, nominally identified by z c , on the figure. The global similarity solution assumes that the 
volume of the gas phase above the interface and below the z = 0 dashed line at t = 0 is 
approximately if not precisely equal to the volume of vapor above the bulk meniscus and below 
the z = Zb line at t = 2.23s, In this way the similarity solution is only concerned with the region Zb 
< z < Zripj where the assumptions of parallel flow and 1-D interface curvature apply. This 
argument serves also to explain the discontinuous nature of the bulk interface at Zb depicted in 
Figure 4. 

For certain container cross-sections of extreme aspect ratio, no value for Tf^. may be computed. 
For example, the 6.5-173.5 equilateral rhombus in Table 1 is the last rhombic section that will 
permit a solution for Tf^ r satisfying Eq. (20). In the limit of 0Cj « 1 and cos# - sinor ; « 1, 
1C*, i becomes increasingly large and Zb moves sufficiently upstream to the point where the 

interfaces in each comer join to form the inscribed circle to the container. At this point the global 
similarity solution is unable to compute Zb (though L(t) and Q{t) remain valid) and Zb and L are of 
equal order. This implies that the bulk meniscus recede rate is comparable to the tip rise rate. It is 
conjectured that such limiting vessels possess a region dominated by curvature associated with 
the inscribed circle with communication between comers, and where the interface may be 
influenced by the Rayleigh instability at long times. No experimental support to this claim is 
offered. 

For the majority of possible containers, solution to the common origin problem provides the 
common value of Zb for all comers. Thus, it is possible to compute the dimensional capillary 
under-pressure at Zb in each comer, namely 

P = 0-345ct ( 09 ) 

J RFKBj) 

where Tf^. is the lower limit of the integral in Eq. (20). Differences in Pj values at Zb between 

comers can provide some insight into the nature of the common origin solution. Obviously, the 
common origin assumption is correct for regular /t-gons and rectangles since the flow in each 
comer is identical. Recalling that for sufficiently negative if, F+ = const., and Pj of Eq. (29) 
approaches a constant at Zb for all comers; a further validation of the common origin assumption 

for containers where I Tf£ B. I » 1 is satisfied in each comer. An example of such a container is 

approached by the geometry of the 6.5-173.5 rhombus in Table 1. However, for other container 
cross-sections, the individual comer capillary under-pressures at Zb differ by as much as the 
theoretical maximum of 26%. If the common origin assumption is correct for these flows too, the 
differences in pressures at Zb between comers must be accommodated within the region of bulk 
meniscus curvature, which has been shrunk to a zero thickness region as sketched in Figure 4 and 
discussed above. 

The analytical approach applies in general if the Concus-Finn condition 0 < nil - otj is satisfied 
in at least one interior comer of the container. The special case of the Concus-Finn condition 
being met in all comers is presented above and supported by experiments performed herein. 
However, it is important to note that for systems exhibiting partial wetting, 9 } > 0, contact angle 
hysteresis may result in significant departures from the above predictions. Also, the analysis at 
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present is limited to polygonal containers without re-entrant comers, though containers with 
re-entrant comers and smoothly curved portions may be treated with a similar approach. 

Drop Tower Experiment s 

Simple experiments to corroborate the theoretical predictions are performed using the 2.2s drop 
tower at NASA’s Glenn Research Center. Partially filled cylinders of triangular cross-section are 
secured to a drop frame, backlit by diffuse light source, and photographed at long working 
distance by a BAV digital video recorder. Release of the experiment package into free fall signals 
the onset of the spontaneous “capillary rise” (imbibition) along the interior comers of the 
container. Sample video images are provided in Figure 7 for a typical drop test, where the 
wetting liquid is observed to redistribute along the comers of the container with the step- 
reduction of gravity, see also Figure 6. The difference in wetting rates due to the different comer 
angles is obvious from such tests. The test fluids used are perfectly wetting liquids on acrylic, 
2 and 5sc Si oil, and the relevant fluid properties are included in Table 2. 

Test Cell and Test Description 

Two test-case prismatic vessels are fabricated for the experiments: 30-60-90 and 45-45-90 right 
triangles. The vessels are fabricated using precision-machined acrylic sections that are carefully 
fused then annealed. To ensure accurate interface height measurements perpendicular to the 
comer axis with minimal though calculable corrections for optical distortion, the test vessels are 
constructed such that the plane bisecting the comer angle of concern is perpendicular to the 
camera line of sight. Figure 8 illustrates the cross-section of the cells tested and provides some 
dimensions. 

Figure 8 also indicates the fused surfaces necessary for vessel fabrication. The 90° comers are 
milled using a zero radius tool. Two each of the vessels are fabricated, the nominal details of 
which are included in Table 2. Some of the experiments are performed at low magnification to 
determine liquid column length (tip length) as a function of time L(/), i.e. Figure 7. Other tests at 
higher magnification are performed to determine the bulk meniscus location Z(t) and transient 
surface elevations h{zJ), i.e. Figure 6. Each test is repeated up to 7 times to verify repeatability as 
well as uniformity of test cell dimensions and comer quality. The Tracker Image Analysis 
System developed by NASA 8 is used to digitize the video images. Where possible an intensity 
threshold algorithm is employed to automatically identify and track the interface frame by frame. 
Measurement accuracy could be determined to be within ±0.3 8mm for the low magnification L{t) 
tests and within ±0.12mm for the higher magnification Z(t) and h(z,t) tests. 

Results and Discussion 


Raw data for the comer tip location L(t) and bulk meniscus location Z(t) are plotted in Figure 9 
for 5cs Si oil in the 30-60-90 right triangular vessel. The initial c-coordinate location of the bulk 
meniscus at t = 0, (see Figure 6), is used as the reference for these measurements. Two tests 
each are provided on the plot to illustrate repeatability. The differences in tip location for the 
different comers of the vessel are plain to see and quantify. 

In Figure 10 through Figure 15, L(t ) or Z(t) is plotted against t U2 as suggested by theory. The 
linear nature of these figures is reassuring. In Figure 10 and Figure 11, L(t) is presented as a 
function of comer angle for the 30-60-90 triangle for 2 and 5cs Si oil, respectively. In Figure 12 
and Figure 13, L(t ) is presented as a function of comer angle for the 45-45-90 triangle for 2 and 
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5cs Si oil, respectively. In all four figures the theoretically determined rise rate K L} from Eq. (23) 
is superposed using a dashed line for comparison to the experiments. The level of agreement 
observed fortifies the functional dependence of the untested parameters captured by the analysis 
leading to eqs. (22) and (23); namely, container size, shape, and surface tension, in addition to 
the tested parameters of comer angle, fluid viscosity, and time. 

Acceptable agreement is achieved in all cases, but appears to worsen in the double limit of 
increased comer angle and time. Experimental uncertainty is suspected for the observed “falling 
off’ of the experimental data from the theoretical slope. For comer angles above 30°, as the 
comer angle increases, the amount of liquid rising in the comer rapidly decreases to the image 
resolution at the tip. Combined with increased optical distortion for the larger angles, the tip 
location measurement technique employed becomes increasingly difficult as the liquid column 
lengthens. This is not the case at short times where there is no ambiguity in tip location, which 
perhaps best explains the agreement of K L with the data at short times for the larger comer 
angles. The fact that the 90° comers are machined rather than fused may also play a role 
retarding the flow as h becomes vanishingly small at the tip. The smaller the comer angle the 
more accurate the data and the better the agreement between K L and the data over the entire test 
duration. Computed and experimentally determined values for K L are listed for comparison in 
Table 4. Nominal values differ by less than ±6%. 

In Figure 14 and Figure 15, Z(t ) is presented vs. t m for the 30-60-90 and 45-45-90 triangles for 
2 and 5cs Si oil. Multiple data sets using shaded and unshaded symbols are presented to indicate 
repeatability. K-j, of Eq. (22) is superposed on the figures using dashed lines and in all cases is in 
favorable agreement with the data as listed in Table 3. Because Zb(t) = Z(0 only for large t, Zb(t ) 
predicted by the global similarity analysis requires longer times before agreement with the data 
than the L(t) predictions. 

The level of agreement shown in Figure 10 through Figure 15 and Table 3 and 4 is a positive 
argument supporting the simplifying theoretical assumption of a common pressure (constant 
height) location for all comers which serves as the origin for all coordinate systems. 

Further data along these lines is presented in Figure 16 and Figure 17 for the 45-45-90 triangle, 
2cs fluid, where h(z,J) is plotted for a variety axial locations z, for both 45° and 90° comers, 
respectively. The z-coordinate location of h(z„t) is referenced to the location of the bulk 
meniscus centerline Zo at / = 0, Points zi above and below the true constant height location are 
selected for digitization. h(Zi,t) values below the constant height location decrease in time to H, 
while h(Zi.t ) values above the constant height location increase in time to H,. The z-coordinate 
location of //, may be determined by interpolation to empirically identify the h{z»t) curve that 
most rapidly achieves a constant value in time. 

Unfortunately, these data suffer increased uncertainty due to pixel resolution limits, especially in 
identifying the location of the comer axis. However, this difficulty introduces an offset to the 
family of curves without affecting the relative relationship of one h(zi,t) curve to another, though 
the slopes at large times are also somewhat difficult to compute. The theoretical constant height 
value computed from Eq. (13) is included as a dashed horizontal line in each figure and, in part, 
indicates the magnitude of the offset which appears small in Figure 16 and Figure 17. 
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Thus, it is possible to determine a range of possible values for the coordinate origin Zh for both 
45° and 90° comers. For example, from Figure 16 it can be discerned that in general 3.50 < zh< 
4.56mm for the 45° comer, whereas from Figure 17, 3.28 < Zh< 4.07mm for the 90° comer. This 
at least implies that differences in Z H values between comers, if any, are small < 0{ 1.3mm) and 
certainly negligible compared to the overall liquid column length > O(50mm). This result also 
supports the use of the significantly simplifying assumption of a common ^-coordinate origin for 
each of the comer flows. 


Concluding Remarks 

The results of an analysis predicting low-g capillary flows in containers with interior comers is 
presented that is based on parallel flow along any interior comer of the container satisfying the 
Concus-Finn wetting condition. The results apply to a large class of cylindrical containers with at 
least one interior comer where the global capillary curvature may be computed via the technique 
of de Lazzer et al. 6 , i.e. cylinders of polygonal section that are devoid of re-entrant comers. The 
strength of the global solution is that the complexities of solving both for the flow and interface 
shape in the region of the bulk interface are avoided due to the highly accurate comer flow 
solutions that control the liquid removal rate from the bulk. 

Closed form expressions are presented to predict important features of the flow such as liquid 
column length in each comer Eq. (23), comer and/or total container flow rate Eq. (12), and 
meniscus recede rate Eq. (22). Fully transient 3-D interfaces may also be quickly computed, see 
eqs. (26)-(28) and Figure 5. Drop tower experiments conducted to verify the theoretical 
predictions vary interior comer angle and fluid viscosity in 30-60-90 and 45-45-90 right 
triangular cylindrical containers. The nominally ±6% agreement between the simple hand- or 
spreadsheet-calculated quantities and the experimental data suggests that the untested functional 
dependence of container size, shape, and surface tension are also correctly captured by the 
analysis, see Table 3 and 4. 

The assumption of a common coordinate origin for each comer flow enables the global similarity 
solution. Comparison with experiments reveals that the resulting flows are rather insensitive to 
this assumption. The closed form expressions derived and benchmarked herein are suitable tools 
for efficient low-g fluids system design such as for the prediction of tank filling characteristics, 
the reorientation and settling transients in vaned containers, and the distribution of flows within 
containers with a variety of interior comers. The results are also applicable to terrestrial 
applications where the influence of gravity is small; i.e. porous wick structures in heat pipes, 
capillary pumped loops. Valuable insights concerning optimum container/pore design may 
be determined in part by inspection of the closed form results. For example, it is shown that a 
30-150 equilateral rhombic section/pore is capable of 9-times the flow capacity of a square 
section of identical area. 
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Figure 4. Notation for global similarity solution in irregular n-gon. 
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Figure 5. Computed surfaces employing eqs. (26) through (28): Left, 30-60-90 and Right, 45-45- 
90 right triangles. Surfaces at time 0, 0.05, 0.3 and 0.6s are shown for 2cs Si Oil in containers 
listed in Table 1. 



Figure 6. Capillary rise in 30° comer of 30-60-90 container, 5cs Si oil ( t = 0, 0.37, 2.23s). Vessel 
is oriented as the vessel in the lower left comer of Figure 8. Notation identified for discussion of 
global similarity solution. 
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Figure 7. Drop tower test: capillary rise in 30-60-90 triangle at times 0, 0.23, 0.7 and 1.63s. 
Quantitative test for 90° comer, right side of image, 5cs Si oil. 


\ 


NAS A/CR— 200 1-21 0900 


16 



Line of sight | 


Figure 8. Test cell sections: quantitative test 
angle noted (dimensions in mm). 
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10. L(t) vs. t 1 ' 2 in 30-60-90 with 2cs. 



Figure 11. L(t) vs. t 112 in 30-60-90 with 5cs. 
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Figure 13. L{t) vs. t m for 45-45-90, 5cs. 
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Table 1. Computed values for Tj^ r for 6=0. Q gC om listed for A fl = 0.5mm* 


Container cross-section Q geom 

K n b + r 

Eq. (25) Eq. (20) Eq.(21) ! 

Triangle 

30-30-120 

1.247 

-0.03498 

-0.03400 

Triangle 

30-60-90 

0.891 

-0.02243 

-0.02202 

Triangle 

30-75-75 

0.862 

-0.02152 

-0.02134 

Triangle 

36.9-53.1-90 

0.742 

-0.01802 

-0.01774 

Triangle 

45-45-90 

0.683 

-0.01638 

-0.01614 

Triangle 

60-60-60 

0.470 

-0.01068 

-0.01056 

Eq. Triangle 


0.470 

-0.01068 

-0.01056 

Square 


0.134 

-0.002741 

-0.002721 

reg-pentagon 


0.0537 

-0.001052 

-0.001045 

reg-hexagon 


0.0256 

-0.000493 

-0.000489 

reg-octagon 


0.00805 

-0.000153 

-0.000152 

Eq. Rhombus 

90-90 

0.134 

-0.002741 

-0.002721 

Eq. Rhombus 

60-120 

0.346 

-0.007495 

-0.007423 

Eq. Rhombus 

45-135 

0.671 

-0.01592 

-0.01569 

Eq. Rhombus 

30-150 

1.249 

-0.03495 

-0.03399 

Eq. Rhombus 

15-165 

2.266 

-0.09413 

-0.08512 

Eq. Rhombus 

6.5-173.5 

3.077 

-0.3642-0.1819 


Eq. Rhombus 

0.1-179.9 

2.26 

— 

— 

Square 


0.134 

-0.002740 

-0.002721 

2:1 Rectangle 


0.114 

-0.002380 

-0.002366 

4:1 Rectangle 


0.0727 

-0.001635 

-0.001623 


Table 2. Nominal test containers and fluid 
properties. D/’s are face widths. 


Cross- 


d 2 

d 3 

Height 

Section 

mm 

mm 

mm 

mm 

30-60-90 

9.8 

17.0 

19.6 

150 

45-45-90 

9.8 

9.8 

13.6 

150 

Test Fluid/ 

V 

<7 

e 

n d 

Vessel matl. 

kg/m s 

N/m 

deg. 

- 

2cs Si oil 

0.00174 

0.0187 

0 

1.390 

5cs Si oil 

0.00465 

0.0197 

0 

1.396 

Acrylic 




1.491 
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Table 3. Measured and predicted bulk meniscus recede 
rates for 30-60-90 and 45-45-90 vessels (mm s 1/L ). 


Vessel 

Fluid 

K z 

exp. 

K z 

Eq. (22) 

30-60-90 

2cs 

-5.47 

- 5.67 


5cs 

- 3.76 

-3.60 

45-45-90 

2cs 

-3.3 

-3.69 


5cs 

-2.1 

-2.34 


Table 4. Measured and predicted liquid 
column rise rates for given comer. 


Vessel/Fluid 

Comer 

k l 

Klj 


angle 

exp. 

theo. 

30-60-90/2cs 

30" 

75.3 

73.3 


60 

50.8 

49.4 


90 

28.6 

28.9 

30-60-90/5cs 

30 

46.4 

46.5 


60 

50.8 

49.4 


90 

28.6 

28.9 

45-45-90/2cs 

45 

47.2 

49.2 


90 

21.5 

23.2 

45-45-90/5cs 

45 

31.6 

31.2 


90 

14.0 

14.8 
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